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Abstract 


A  forced-vibration,  non-resonant  method  for  the  determination  of  the  dynamic 
mechanical  properties  of  elastomeric  materials  in  the  audio  frequency  range  was  developed  and 
evaluated  using  a  commercial  damping  material.  Isodamp  €-1002  from  EAR  Specialty 
Composities.  Measurements  of  storage  modulus  and  loss  factor  were  made  in  the  1 1  Hz  to  10  kHz 
range  on  the  damping  material  using  compressive  and  tensile  excitation.  The  elastomer  had  a  loss 
factor  near  1.0  and  a  storage  modulus  which  varied  between  6  x  10^  Pa  and  3  x  10*  Pa  in  the 
frequency  range  studied.  For  tensile  excitation,  the  error  in  the  calculated  modulus  was  low  up  to 
the  resonance  frequency  of  the  sample,  usually  ~1700  Hz.  For  compressive  excitation  of  thin 
sheets,  the  error  was  greatest  at  the  low  and  high  ends  of  the  frequency  range. 


Resume 


Une  m6thode  exploitant  des  vibrations  forc^es  non  rdsonantes  pour  ddterminer  les 
propri^tes  mechaniques  dynamiques  des  elastomeres  dans  la  gamme  des  frequences  audibles  a  6t6 
mise  au  point  et  6valuee  au  moyen  d'un  mat6riau  amortisseur  commercial,  I'lsodamp  C-1002 
d'EAR  Specialty  Composites.  Le  module  de  conservation  et  le  facteur  de  pertes  ont  6\6  mesurds, 
sur  I'intervalle  de  1 1  Hz  a  10  kHz,  dans  le  mat6teriau  amortisseur  soumis  ^  des  contraintes  tantot 
de  compression,  tantot  de  traction.  L'dastomere  a  pr^sent^  un  facteur  de  pertes  de  pres  de  1, 0  et 
un  module  de  conservation  variant  entre  6  x  10^  et  3  x  10^  Pa  dans  la  gamme  de  frequence  donne. 
En  Traction,  I'erreur  sur  le  module  calculi  a  6t6  faible  jusqu'a  la  frequence  de  resonance  de 
I'eprouvette,  normalement  ~1700  Hz  environ.  En  compression,  I'erreur  pour  des  feuilles  minces  a 
6te  la  plus  61ev6e  aux  valeurs  maximale  et  minimale  de  la  gamme  de  frequences. 
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Introduction 


The  unique  properties  of  polymers  make  them  ideally  suited  for  damping  structure- 
home  vibrations  when  used  in  anti-vibration  mounts  or  sheet  damping  applications,  and  in 
reducing  air  or  water-borne  sound  when  used  in  anechoic  coatings^  The  frequency  and 
temperature  dependence  of  the  viscoelastic  properties  of  polymers  is  essential  information 
for  the  selection  and  design  of  materials  for  noise  and  vibration  damping.  There  are  a 
variety  of  commercially  available  instruments,  known  as  Dynamic  Mechanical  Thermal 
Analyzers  (DMTA),  that  are  able  to  measure  the  complex  modulus  of  polymers  over  a  wide 
temperature  range  (-100  to  5(X)  oC);  however,  the  frequency  range  for  these  instruments 
is  limited  to  lower  frequencies,  typically  lO'^  Hz  to  10  Hz.  If  data  at  higher  frequencies  are 
required,  then  an  apparatus  is  usually  constructed,  consisting  of  a  rig  to  clamp  the  sample, 
transducers  to  measure  force  and  displacement,  and  electronics  to  process  the  signals. 

An  important  consideration  in  the  design  of  a  complex  modulus  apparatus  is  the 
excitation  technique.  Resonance  methods,  in  which  at  least  one  end  of  a  sample  can  freely 
move,  are  not  well  suited  to  high  loss,  low  stiffness  samples^  and  data  is  usually  analyzed 
at  only  a  few  frequencies  corresponding  to  the  resonance  frequencies  of  the  sample.  Forced 
vibration  non-resonant  excitation,  in  which  the  sample  is  clamped  at  one  end  and  driven  at 
the  other,  is  well  suited  to  low  stiffness  samples  and  yields  data  over  a  continuous  range  of 
frequencies  for  a  given  sample. 

Another  important  aspect  of  the  excitation  technique  is  the  mode  of  excitation. 
Samples  may  be  excited  in  compression,  tension,  torsion,  shear,  or  flexure.  The  choice  of 
excitation  mode  depends  on  the  stiffness  of  the  sample,  and  whether  shear.  Young's 
(extensional),  longitudinal,  or  bulk  modulus  is  of  most  interest.  For  elastomers  in  their 
rubbery  state,  shear  or  tensile  excitation  works  well,  but  stiffer  samples  may  require  a 
bending  mode.  A  major  disadvantage  in  a  bending  experiment  is  that  the  Poisson's  ratio  of 
the  polymer,  which  can  vary  from  0.3  to  0.5,  should  be  known  in  order  to  accurately 
calculate  shear  and  Young's  moduli. 

In  this  paper,  the  forced-vibration  non-resonant  excitation  method  is  described  in 
detail  for  two  experimental  configurations. 

Experimental  Procedure 

Measurements  were  carried  out  on  a  commercial  damping  polymer.  Isodamp  C- 
1(X)2  from  EAR  Specialty  Composites.  Soxhlet  extraction  of  this  polymer  with  acetone 
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showed  it  to  contain  about  50  wt%  liquid  plasticizer.  Pyrolysis/Gas  Chromatography- 
Mass  Spectrometry  (GC-MS)  and  Infrared  (IR)  analysis  identified  the  plasticizer  as 
dioctylphthalate  and  the  deplasdcized  polymer  as  polyvinyl  chloride  (PVC),  or  a  PVC 
copolymer.  Powder  x-ray  diffraction  and  Energy  Dispersive  X-ray  Analysis  (EDXA)  of 
the  polymer  revealed  that  the  polymer  contained  Ti02  as  an  inorganic  filler. 

The  complex  Young's  modulus  of  the  EAR  material  was  measured  using  two 
ctmfigurations: 

Configuration  I  Figures  1-4  illustrate  the  experimental  set  up  for  tensile  excitation  of  a 
rectangular  polymer  sample  (3.08  mm  x  11.80  mm  x  26.61  mm).  The  force  transducer  (B 
&  K  82(X))  was  bolted  into  a  massive  copper  plate  (5  cm  x  38  cm  x  71  cm)  that  served  as  a 
rigid  support  for  this  system.  The  top  and  bottom  of  the  sample  were  held  in  place  by 
hand-tightened  clamps.  The  top  clamp  remained  stationary,  while  the  bottom  clamp  was 
driven  with  broadband  excitation  provided  by  a  B  &  K  4809  Vibration  Exciter,  which  was 
supported  by  an  adjustable  lab  jack.  The  motion  of  the  driven  end  of  the  sample  was 
monitored  by  an  accelerometer  (B  &  K  4374).  The  force  (F)  and  acceleration  (x)  signals 
were  amplified  by  charge  amplifiers  (B  &  K  2635)  before  being  frequency  analyzed  by  a 
digital  filter  analyzer  (B  &  K  2133)  in  1/12  octave  mode.  Both  autospectra  and  cross 

spectra  were  collected  to  obtain  the  magnitudes  k  and  F,  as  well  as  the  phase  angle  (9) 
between  them.  The  spectral  data  were  transferred  over  an  IEEE  bus  to  a  Jemini  286 
computer.  Calculation  of  storage  moduli  and  loss  factor  were  carried  out  using  Equations  6 
and  8  (see  Results  and  Discussion). 

Configuration  II  Figures  5  and  6  illustrate  the  experimental  set  up  for  compressive 
excitation  of  a  thin  sheet  (3.03  mm).  The  force  transducer  was  bolted  directly  to  the 
shaker  and  to  a  15  mm  x  16  mm  aluminum  block  that  served  to  transmit  the  force  from  the 
vibration  exciter  to  the  sample,  which  was  compressed  between  the  aluminum  block  and 
the  rigid  support.  Two  accelerometers  were  used  to  monitor  the  motion  of  the  vibration 
exciter  (and  aluminum  block)  and  the  motion  of  the  nominally  rigid  support.  As  shown  in 
Bgure  6,  the  amplified  accelerometer  signals  were  differenced  by  a  differential  ac  amplifier 
(Brookdeal  9454).  The  net  acceleration  (j)  and  force  (F)  signals  were  then  frequency 
analyzed  in  1/12  octave  mode.  Both  autospectra  and  cross  spectra  were  collected  to  obtain 

the  magnitudes  y  and  F,  as  well  as  the  phase  angle  (9)  between  them.  Calculation  of 
storage  moduli  and  loss  factor  were  carried  out  using  Equations  18  and  20  (see  Results  and 
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Discussion).  The  spectra  for  this  conflguration  were  smoothed  using  a  second  order  nine 
point  Savitzky-Go  ay  algorithm^.  There  was  no  need  to  smooth  data  from  Conflguradon  I. 

Results  and  DiLcussion 


Configuration  I  Fi||ures  7  and  8  show  the  measured  force  (F),  acceleration  (^,  and  phase 
angle  (q>)  for  an  elastomer  sample  excited  in  tension.  Notice  that  the  phase  angle  is  a 


smooth  function  untir,  {le  first  resonance  peak  occurs  at  ~  1700  Hz.  In  contrast,  the  phase 
angle  for  Configuratio^JI  was  well  behaved  at  high  frequency,  but  poorly  behaved  at 
frequencies  below  100 1%  For  Configuration  I,  the  force  transducer  was  stationary  and 
the  only  force  acting  on  the  transducer  was  from  tensile  excitation  of  the  sample.  Since  the 
strain  on  the  sample  was  smaW'^  (<  0.4  %),  the  force  (F)  is  proportional  to  the  displacement 
(x)  as  described  by  Hooke's  Law:  \ 

F  =  kx  [1] 

where  k  is  the  sample  stiffness.  F  and  x  are  time  (t)  varying  quantities  and  thus  have 
amplitude  and  phase: 

x  =  xexp[ia)t]  [2] 

F  =  Fexp[i((ot  +  <p)]  [3] 

In  the  above  equation,  the  barred  symbols  represent  amplitudes  of  complex  variables,  co  = 
2jcf  is  the  angular  frequency,  <p  is  the  measured  phase  angle  between  F  and  x,  and  i  =  . 

Substituting  [2]  and  [3]  into  [1]  gives 

k  =  (F/ x)  exp[i<p]  [4] 

For  a  tensile  experiment  on  a  sample  of  length  L  and  cross  sectional  area  A,  the  Young's 
modulus  is  defined^  to  be  E  =  k(L/A).  Noting  from  [2]  that  x  =  -co^x,  it  follows  that 


E  =  -  t  I  co2  exp[i<p]  [5] 

The  real  and  imaginary  parts  of  the  complex  modulus,  or  storage  and  loss  modulus,  are 


E  =  -L 


— f  ©2  cos  9 
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m 


E'’  = 


L  J  FI  . 

X  ”  ~^{  sin  (p 

A  [  X  J 


The  material  loss  factor,  11,  is  related  to  the  phase  angle  between  stress  and  strain  on  the 
sample,  S: 

E'* 

Ti  =  tan  5  =  =  tan  (p  [8] 

In  this  case  the  measured  phase  angle,  (p,  is  identical  with  5,  the  angle  between  stress  and 
strain. 

The  storage  moduli  and  loss  factors  (Figures  9 A  and  lOA)  were  calculated  by 
Equations  6  and  8,  using  data  shown  in  Figures  7A,  7C,  and  8A.  In  this  configuration,  the 
resonance  at  17(X)  Hz  results  in  large  errors  in  the  calculated  storage  moduli  and  loss 
factors  at  frequencies  approaching  resonance  and  above  resonance.  The  modulus  increases 
monotonically  with  frequency,  with  rubber-like  values  (<10^  Pa)  at  frequencies  up  to  100 
Hz,  and  with  values  characteristic  of  the  glass  transition  (-*  10^  Pa)  at  higher  frequencies. 
The  calculated  loss  factors  as  a  function  of  frequency  forms  a  smooth  curve  up  to  300  Hz, 
after  which  the  resonance  causes  large  errors. 

Curve  B  in  Figures  9  and  10  represents  data  from  a  previous  study  on  Isodamp  C- 
1002  by  Sun  and  Mitchell^.  The  technique  used  by  Sun  and  Mitchell  is  a  transfer  function 
method  with  an  associated  "lumped  mass"  model.  Calculation  of  moduli  with  the  latter 
method  requires  the  solution  of  several  non-linear  equations  for  each  frequency  of  interest, 
in  contrast  to  the  method  described  here,  for  which  post-processing  of  experimental  data  is 
quite  straightforward.  Agreement  between  the  two  techniques  is  reasonable  up  to  the 
resonance  frequency  of  the  sample.  Adjustment  of  pre-load  to  near  zero  values  was  carried 
out  by  adjusting  the  height  of  the  vibration  exciter  with  the  lab  jack.  Increasing  the  pre-load 
resulted  in  larger  calculated  moduli,  but  a  more  quantitative  study  of  this  relationship  will 
be  reserved  for  future  experiments  in  which  a  static  load  cell  will  be  incorporated  into  the 
system. 

Configuration  11  For  this  set  up,  the  use  of  a  second  accelerometer  was  found  to  give 

a  smoother  phase  angle  curve.  Figures  7  and  8  show  the  measured  force  (F),  net 
acceleration  ^),  and  phase  angle  (tp)  for  an  elastomer  sheet  excited  in  compression  against 
a  rigid  support.  This  type  of  excitation  required  larger  forces  for  smaller  displacements 
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than  Configuration  I,  since  the  stress  was  across  the  sample  thickness,  not  length.  The  low 
signal/noise  ratio  may  have  contributed  to  errors  in  the  phase  angle  at  frequencies  below 
100  Hz,  which  resulted  in  errors  in  the  calculated  moduli  (see  below). 

The  force  acting  on  the  transducer  included  the  inertial  force  required  to  move  the 
aluminum  block.  Thus,  the  equation  of  motion  for  this  system  must  include  a  mass  term^: 

F  =  mx  +  k(H-Xo)  [9] 

where  F  is  the  measured  force,  m  is  the  mass  of  the  sample  block,  x  is  the  displacement  of 
the  sample  block,  Xq  is  the  displacement  of  the  support,  and  k  is  the  stiffness  of  the 
sample.  With  the  exception  of  mass,  all  variables  in  [9]  are  time-varying  quantities  and  thus 
have  an  amplitude  and  phase.  If  we  let  y  =  x  -  Xq  be  the  net  displacement,  then  y,  F,  and 
Xo  may  be  expressed  as 

y  =  yexp[icDt]  [10] 

F  =  F  exp[i(cot  +  <p)]  [11] 

Xo  =  Xo  exp[i(cot  +  0)]  [1 2] 

Substituting  y  =  x  -  Xo  and  k  =  y  +  Xo  into  [9]  gives 

F=my+mxo  +  ky  [13] 

Noting  from  [10]  that  y  =  -CO^y, 

F=my+mxo--^k  [14] 

0)2 

k=  -  -  0)2  +  o)2mf  1  +  [15] 

y  {  y ) 

If  the  support  is  very  rigid,  then  Xq  «  y  and  Equation  [15]  can  be  approximated^  by 
k=  -?o)2+  o)2m  [16] 

y 

The  validity  of  this  approximation  requires  verification  by  observation  of  the 
accelerometer  signals.  Substituting  [10]  and  [11]  into  [16]  gives 
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[17] 


k=  — ^  (1)2  exp[i(p]  +  0)2m 

y 

For  a  sample  of  thickness  d  and  cross  sectional  area  A,  the  extensional  modulus  E  =  kd/A. 
It  follows  from  [17]  that  the  real  and  imaginary  parts  of  the  complex  modulus  are 


E*  = 


d 

A 


0)2  COS  <p  +  (o2m 


[18] 


E"  = 


d 

A 


0)2  sin  <p 


[19] 


In  this  case,  the  measured  phase  angle,  (p,  is  different  from  the  angle  between  stress  and 
strain  on  the  sample,  5,  and  so  tan  5  ^  tan  9: 


F  5  . 

-  -3—  (1)^  sin  (p 

E"  V 

“n  =  tan  5  =  ^ - 

F 

-  -::r-  (o2  cos  9  +  a)2m 

y 


[20] 


E" 

At  low  frequencies.  Equation  [20]  reduces  to  -gr  =  tan  9. 

In  order  to  calculate  the  inertial  mass  of  the  system  to  be  used  in  Equations  18  and 
20,  the  aluminum  block  is  excited  without  the  sample  in  a  free  vibration  experiment.  The 
effective  mass  is  simply  m  =  F/3c,  and  as  shown  in  Figure  1 1,  is  independent  of  frequency. 

Figure  12,  curve  A,  shows  the  storage  modulus  as  a  function  of  frequency  as 
determined  using  Configuration  II  and  Equation  18  for  Isodamp  C-1002.  There  is  excellent 
agreement  with  the  literature  data  (curve  B)  in  the  range  300  Hz  -  60(X)  Hz.  The  loss  factors 
calculated  from  Equation  20  do  not  compare  as  well  to  the  literature  data  as  the  storage 
moduli  do,  and  even  after  smoothing  the  loss  factor  spectrum  appears  noisy  (Figure  13). 
The  errors  associated  with  this  method  seem  to  be  related  to  errors  in  the  phase  angle  at  low 
frequencies  and  to  the  0)2m  term  at  high  frequencies.  This  method  was  particularly 
sensitive  to  the  amount  of  static  pre-load,  with  the  most  consistent  results  obtained  at  near 
zero  pre-loads.  Adjusting  the  load  was  sometimes  difficult,  since  the  sample  was  held  in 
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place  by  compressing  it  between  the  aluminum  block  and  the  support,  and  this  necessitated 
a  pre-strain  greater  than  the  dynamic  strain.  An  increase  in  the  static  pre-load  generally 
resulted  in  increased  calculated  moduli. 

The  experimental  set  up  for  compression  may  be  improved  by  using  an  impedance 
head  in  place  of  separate  transducers,  and  a  mass  compensation  circuit  to  more  accurately 
account  for  the  mass.  The  impedance  head  combines  a  force  transducer  and  accelerometer 
in  a  single  housing,  allowing  the  measurement  of  impedance  at  the  driving  point  of  the 
transducer. 

Configuration  I  was  found  to  be  a  more  reproducible  technique  than 
Configuration  II,  with  better  signal/noise  ratio  and  less  sensitivity  to  sample  positioning. 
Also,  adjustment  of  pre-load  was  more  straightforward  for  Configuration  I.  The  two 
experimental  configurations  described  here  are  complementary  in  that  they  are  applicable  to 
different  frequency  ranges. 

Conclusion 

The  measurement  of  complex  Young's  modulus  over  a  range  of  audio  frequencies 
was  demonstrated  on  a  commercial  damping  material  using  two  variations  of  a  forced 
vibration,  non-resonant  technique.  For  tensile  excitation,  a  simple  model  was  successful  in 
interpreting  the  data  up  to  the  resonance  frequency  of  the  sample.  Compressive  excitation 
required  a  model  that  incorporated  an  inertial  mass  acting  on  the  force  transducer,  but  there 
was  no  interference  from  sample  resonance.  Sources  of  error  in  both  ca^es  could  be  traced 
to  the  phase  angle,  which  is  sensitive  to  resonance  phenomena  and  low  signal  to  noise 
ratio. 

The  method  described  here  has  several  advantages  over  other  techniques  for 
measuring  complex  modulus  of  elastomers.  First,  because  it  is  a  non-resonant  method  and 
broadband  excitation  is  used  in  conjunction  with  spectral  analysis,  data  is  available  over  a 
large  range  of  frequencies,  whereas  resonance  methods  are  limited  to  just  a  few 
frequencies.  Second,  the  stiffness  of  the  sample  is  directly  measured,  which  simplifies 
data  analysis  and  reduces  calculation  errors  associated  with  methods  that  utilize  transfer 
impedance  measurements.  Finally,  unlike  the  double-lap  shear  method  or  the  transfer 
function  method,  no  adhesive  was  used  in  clamping  the  sample,  which  is  advantageous  if 
large  numbers  of  samples  are  to  be  tested. 
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RIGID  SUPPORT 


CONFIGURATION  I 


Figure  1.  Schematic  diagram  of  set  up  for  tensile  excitation  of  polymer  samples. 
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Figure  2.  Photograph  of  experimental  apparatus.  (A)  plotter,  (B)  dual  channel  spectrum 
analyzer,  (C)  power  amplifier  for  vibration  exciter,  (D)  charge  amplifiers,  (E)  copper  plate, 
(F)  vibration  exciter.  The  vibration  exciter  is  supported  by  an  adjustable  lab  jack. 


- 11  - 


Figure  3.  Close  up  of  clamping  arrangement  for  Configuration  I.  Under  label  E;  force 
transducer,  top  clamp,  sample,  bottom  clamp,  accelerometer. 
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Figure  4.  Block  diagram  of  instruments  for  Configuration  I. 


- 13- 


ACCELEROMETER 


RIGID  SUPPORT 


CONFIGURATION  II 


Figure  5.  Schematic  diagram  of  set  up  for  compressive  excitation  of  polymer  sheets. 
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Figure  6.  Block  diagram  of  instruments  for  Configuration  II. 
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Figure  8.  Phase  angle  between  acceleration  and  force  signals  for  EAR  C-1002,  obtained 
from  measurement  of  cross  spectra.  (A)  Configuration  I  and  (B)  Configuration  II. 
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part  of  cxtensional  modulus  versus  frequency  for  EAR  0-1002  (A) 
nd  (B)  From  Reference  5. 
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Figure  10.  Loss  factor  versus  frequency  for  EAR  C-1002  (A)  Configuration  I  and 
(B)  From  Reference  5. 
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Figure  1 1 .  Plot  of  measured  acceleration  and  force,  and  calculated  mass  in  a  free  vibration 
experiment  for  estimation  of  inertial  mass  term  to  be  used  in  Equations  19-20. 
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Figure  13.  Loss  factor  versus  frequency  for  EAR  C-1002  (A)  Configuration  II  and 
(B)  From  Reference  5. 
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